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SECTION  1 
INTRODUCTION 


This  report  describes  work  on  the  WEDCOM  environmental,  propagation  and 
system  models  performed  during  1982.  This  effort  is  a  continuation  of  work  to 
prepare  a  revised  version  of  WEDCOM,  called  WEDCOM  V,  by  January  1984.  The 
environmental  and  system  modeling  tasks  are  new.  The  propagation  model  tasks  were 
concerned  with  the  implementation  and  verification  of  the  models  developed  during 
the  previous  effort,  reported  in  Reference  1. 

Section  2  describes  work  on  the  environmental  models  to  incorporate  current 
models  and  provide  for  code  selection  of  calculation  times  and  location  of  verti¬ 
cal  ionization  profiles  along  the  propagation  path. 

Section  3  describes  the  Implementation  and  modifications  to  the  new  VLF/LF 
propagation  models  described  in  Reference  1.  The  code  structure  was  optimized  to 
allow  efficient  computations  consistent  with  the  propagation  sensitivities  with 
the  objective  to  reduce  the  link  evaluation  time  required  by  WEDCOM  IV.  Signifi¬ 
cant  computational  time  reductions  were  obtained,  but  they  were  case  dependent. 
Computation  time  reductions  achieved  are  typically  factors  of  three  to  five  for 
the  VLF  and  LF  models  and  five  to  ten  for  the  ELF  model.  The  actual  time  reduc¬ 
tion  is  complicated  by  a  new  procedure  which  more  accurately  accounts  for  iono¬ 
spheric  and  ground  conductivity  variations  than  did  the  WEDCOM  IV  code.  The 
number  of  detailed  reflection  coefficient  and  diffraction  computations  may  be 
greater  than  might  have  been  specified  by  the  WEDCOM  IV  user.  However,  the 
overall  code  efficiency  and  propagation  parameter  sensitivity  has  been  improved. 

The  propagation  model  verification  was  accomplished  by  combining  these  models 
with  the  WEDCOM  IV  (Reference  2)  structure  and  environmental  models.  A  determina¬ 
tion  of  model  accuracy  and  computational  time  savings  was  made  by  comparisons  with 
the  WEDCOM  IV  code. 

The  new  ELF  propagation  model  implementation  is  described  in  Section  4. 
Comparisons  to  other  models  are  made  for  model  verification  and  the  new  subrou¬ 
tines  are  described. 


VLF/LF  system  modem  models  developed  by  the  Stanford  Research  Institute 
(References  3  and  4),  were  adapted  for  WEDCOM.  This  effort  is  described  in 
Section  5.  Also  described  is  the  work  on  an  atmospheric  noise  model  that  includes 
the  effects  of  a  nuclear  environment. 


SECTION  2 

ENVIRONMENTAL  MODEL  ADDITIONS 


This  section  describes  the  following  modifications  to  the  WEDCOM  code  en¬ 
vironmental  models: 

1.  Implement  current  environmental  modeling. 

2.  Provide  for  code-selected  computational  times. 

3.  Provide  for  code-selected  vertical  path  locations. 

4.  Provide  for  ionization  file  retention  for  subsequent  system 
computations. 

IMPLEMENT  CURRENT  ENVIRONMENT  MODELS 

The  environment  models  used  in  WEDCOM  IV  are  those  developed  for  WEPH  IV 
Mod  2  in  1976  (Reference  5).  The  following  revised  environment  models  prepared 
for  the  WEPH  code  (Reference  6)  and  the  WESCOM  code  (Reference  7)  have  been 
implemented  into  WEDCOM. 

Model 

Ambient  Atmosphere 
Ambient  Ionosphere 

Atmospheric  Chemistry 
Post-Stabilization  Debris 

High-Altitude  Fireball  Plume 

CODE-SELECTED  CALCULATION  TIMES 

The  WEDCOM  IV  user  must  specify  the  calculation  times  to  be  used  for  system 
evaluation.  This  is  satisfactory  for  users  who  have  some  familiarity  with  the 
intrinsic  characteristics  of  nuclear  environments  and  system  responses.  However, 
adequate  time  selection  can  be  difficult,  particularly  for  complex  scenarios. 


Comments 

ROSCOE  Optics  Model 

Provides  for  variation  with  solar 

and  geographic  conditions 

ROSCOE  Optics  Model  1978 

Improved  agreement  between  analytic 

and  numerical  models 

New  late-time  fireball  and  debris 

models 


In  WEDCOM  V  provision  is  being  made  for  both  user-  and  code-selected  cal¬ 
culation  times.  For  code  selection  the  user  specifies  the  time  interval  of 
interest  and  the  code  will  select  a  set  of  evaluation  times  that  take  into  account 
spatial  and  temporal  considerations. 

A  simplified  flow  chart  showing  the  times  selected  for  phenomenology  calcu¬ 
lations  are  shown  in  Figure  1.  The  output  is  a  phenomenological  file  for  all 
evaluation  times  determined  by  phenomenology  variations.  This  file  provides  the 
burst  characteristics  (location,  regions,  etc). 

A  time  array  (1,  30,  120,  300,  900,  1800,  3600  seconds  and  every  hour  there¬ 
after)  is  defined  beginning  at  each  burst  time.  The  phenomenology  computations 
require  little  execution  time  when  compared  to  that  required  by  the  ionization  and 
propagation  models.  Thus,  for  a  single  burst  the  nuclear  phenomenology  character¬ 
istics  are  determined  after  every  detonation  at  the  above  time  intervals  (within 
the  time  frame  interval).  For  a  multiburst  case  the  evaluation  times  occur  at  the 
same  time  intervals  after  each  subsequent  burst,  but  do  not  include  the  evaluation 
times  associated  with  prior  bursts.  The  phenomenology  characteristics  are  stored 
on  file  for  use  in  determining  vertical  ionization  profiles  at  selected  locations 
between  transmitter  and  receiver.  To  minimize  the  number  of  ionization  and 
propagation  calculations,  calculation  times  for  a  burst  that  is  determined  to  be 
unimportant  (described  below)  are  replaced  with  times  selected  for  prior  important 
bursts  before  determining  the  ionization  environment. 

CODE-SELECTED  VERTICAL  PATH  LOCATIONS 

In  WEDCOM  IV  the  ionization  environment  is  computed  for  verticals  equally 
spaced  along  the  great  circle  path  between  transmitter  and  receiver.  In  many 
cases  this  procedure  selects  too  many  verticals  in  undisturbed  regions  and  too  few 
where  nuclear  disturbances  are  important.  A  new  procedure  has  been  developed  that 
provides  for  code-selected  vertical  path  locations  that  are  determined  on  the 
basis  of  the  disturbed  environment. 

The  procedure  is  based  on  first  determining  the  importance  of  the  prompt  and 
delayed  radiation  from  each  burst  at  vertical  path  locations. 

A  maximum  range  of  effects,  ,  for  every  prompt  or  debris  region  is 

defined.  A  prompt  region  is  not  considered  if  the  detonation  altitude  is 

less  than  25  km  because  the  region  is  too  small  for  propagation  sensitivity 
studies. 


Figure  1.  CONWED  modifications  illustrating  the  logic  for 
selecting  the  initial  computation  times. 


For  higher  detonation  altitudes,  R^  for  prompt  radiation  is  defined  as 


R 


d 


500  km,  Hb  <  85  km 

500  +  R  cos  ([R  +85] /[R  +HJ)  km,  H,  >  85  km 
e  e  e  b  b 


(1) 


where  R  is  the  earth's  radius, 
e 

A  debris  region  is  not  considered  if  the  debris  altitude,  ,  is  below 

15  km.  This  elimination  is  also  due  to  negligible  propagation  sensitivity.  For 
higher  debris  altitudes,  R^  is  defined  as 

!R, .  +  25  H.  -  375  km,  for  H,  <  25  km 

Ida  d 

(2) 

R,  ,  acos  ([R  +25] / [ R  +H,])  +  R.  .  +  250  km,  for  H,  >  25  km  . 

id  e  e  d  Id  d 


where  R,  ,  is  the  radius  of  the  debris. 

Id 

A  simplified  flow  chart  for  the  vertical  path  selection  is  shown  in  Figure  2. 
Beginning  at  the  transmitter,  the  propagation  path  is  incremented  toward  the 
receiver  (nominally  in  500-km  steps) .  At  each  segment  location  the  range  to  the 
prompt  source  or  debris  center  is  determined.  If  this  range  is  smaller  than 
R^  ,  the  segment  and  phenomenological  time  are  flagged.  A  new  phenomenology 
file  is  written  (for  up  to  twenty  flagged  evaluation  times)  which  will  sub¬ 
sequently  be  utilized  for  ionization  calculations. 

Once  the  segments  have  been  flagged,  logical  tests  are  made  to  determine  if 
at  least  two  adjacent  segments  are  flagged  at  a  given  evaluation  time.  If  not, 
the  single  segment  is  ignored  because  the  affected  region  is  small  and  assumed  to 
be  negligible.  However,  if  there  are  two  or  more  adjacent  flagged  segments,  the 
neighboring  segments  are  flagged.  This  will  help  refine  the  disturbed  region 
definition  along  the  propagation  path. 

RETENTION  OF  IONIZATION  FILE 

If  the  time  frame  specified  by  the  user  is  large  compared  to  the  time  for 
signficiant  environmental  variations,  the  number  of  evaluation  times  selected  by 
the  code  may  be  excessive,  in  turn  requiring  long  and  unacceptable  computation 
times.  The  ionization  file  prepared  by  a  WEDCOM  execution  will  be  retained  to  be 
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Figure  2.  Logical  flow  implemented  in  Subroutine  ENVIRN  to  determine 
computation  locations. 


utilized  by  a  subsequent  execution.  Reference  1  suggests  that  this  file  be  used 
to  determine  the  ionization  (by  interpolation)  at  the  user  specified  times.  This 
feature,  while  straightforward  to  implement,  would  require  significant  new  coding. 
Thus  we  have  elected  not  to  include  it  and  to  require  the  user  to  select  from  the 
preselected  computation  times. 

The  ionization  file  might  have  a  large  number  of  evaluation  times,  and  while 
these  times  may  represent  the  environment  sensitivities  they  do  not  necessarily 
correspond  to  the  propagation  or  system  sensitivities.  It  is  planned  to  require 
that  the  initial  execution  provide  output  such  as  ionospheric  absorption,  re¬ 
flection  altitudes,  and  gradients  for  the  systems  being  analyzed.  These  outputs 
require  minimal  computation  resources  and  would  assist  the  user  in  selecting 
evaluation  times  for  a  subsequent  execution. 


SECTION  3 

VLF/LF  PROPAGATION  MODEL  IMPROVEMENTS 


The  effort  reported  here  consists  of  the  modification,  implementation,  and 
testing  of  the  model  techniques  described  in  Reference  1,  Changes  and  improve¬ 
ments  were  made  in  the  overall  structure  of  VLF  and  LF  computations  and  in  various 
submodels  in  both  the  LF  and  VLF  model.  Some  of  the  modeling  changes  apply  to 
both  VLF  and  LF  models.  The  following  subsections  will  parallel  the  topics 
described  in  Sections  2,  3,  and  4,  Reference  1. 

LF  MODEL  FORMULATION 

This  formulation  was  implemented  as  described  in  Subsection  2.2,  Reference  1. 
Figure  3  illustrates  the  logical  flow  in  Subroutine  BHRHS  which  defines  the 
ionospheric  profiles  at  points  along  the  propagation  path.  The  profiles  differ 
from  one  another  by  the  ionization  value  and/or  the  magnetic  propagation  azimuth 
and  dip  angle.  The  profiles  apply  to  path  segments  approximately  equal  to  the 
spacing  between  profiles. 

Reflection  altitudes  for  the  segments  are  not  linearly  interpolated  between 
adjacent  profiles  as  suggested  in  Reference  1,  and  as  was  done  in  past  WEDCOM 
versions.  Studies  of  a  variety  of  reflection  regions  suggested  that  the  inter¬ 
polated  values  did  not  necessarily  provide  more  accurate  modeling.  Therefore,  a 
uniform  reflection  altitude  corresponding  to  the  profile  at  the  segment  center  is 
used  for  each  segment.  This  minimizes  the  number  of  reflection  coefficient 
computations  and  reduces  the  amount  of  programming. 

The  complex-valued  coefficients  for  the  analytic  ionospheric  reflection 
coefficient  formulation  are  computed  as  described  in  Section  3,  Reference  1  and 
illustrated  in  Figure  4.  Note  that  a  maximum  of  three  calls  to  the  reflection 
coefficient  model  are  made.  The  incident  angle  difference  between  the  most 
oblique  and  steepest  Incident  angles  for  each  segment  must  be  greater  than  0.5 
degrees  to  require  two  computations,  and  greater  than  3  degrees  to  require  three 
computations.  Otherwise  a  single  computation  is  performed. 


Figure  3.  Logical  flow  for  defining  ionospheric  segments  along 
the  propagation  path. 


ogical  flow  for  computing  the  complex-valued  A  and  B  coefficients  for 
he  analytic  reflection  coefficient  formulation  used  in  the  LF  model. 
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Sky  Wave  Calculations 

Section  2.4,  Reference  1,  described  the  algorithms  for  the  sky  wave  calcu¬ 
lations.  'The  earth-detached  (Whispering  Gallery)  ray  paths  were  also  described. 
These  are  ray  paths  that  propagate  from  one  ionospheric  reflection  point  to 
another  without  an  earth  reflection,  and  they  can  be  very  important  for  highly 
elevated  antennas.  However,  the  application  of  WEDCOM  seldom  requires  antenna 
altitudes  above  that  of  normal  aircraft;  thus  earth-detached  ray-path  model  is  not 
required  and  was  not  included.  This  omission  Insures  computation  time  require¬ 
ments  consistent  with  the  objectives  of  WEDCOM.  There  are  special  considerations 
for  highly  elevated  antennas  that  are  not  easily  incorporated  into  WEDCOM,  espe¬ 
cially  for  disturbed  environments.  Thus,  WEDCOM  V  will  simply  output  a  flag 
noting  that  a  user  has  specified  an  antenna  altitude  greater  than  10  km. 

Ground  Wave/Earth  Surface  Diffraction 

Sections  2.5  and  2.6,  Reference  1,  described  the  new  ground  wave  and  earth 
surface  diffraction  models  for  WEDCOM  V  implementation.  These  models  are  major 
users  of  computation  time,  second  only  to  the  ionospheric  reflection  coefficient 
model.  Repetitive  computations  were  eliminated  from  the  model,  as  reported  in 
Reference  1.  Further  reductions  in  repetitive  computations  were  possible  by 
utilizing  the  model  characteristic  that  a  maximum  of  six  possible  conductivities 
can  be  output  from  the  worldwide  ground  conductivity  map  (see  Section  2.3,  Refer¬ 
ence  1).  Thus,  the  WEDCOM  V  code  initializes  the  ground  wave/diffraction  models  a 
maximum  of  six  times,  depending  on  the  conductivity  of  segments  of  the  propagation 
path. 

Reflection  Altitude 

An  analytic  reflection  altitude  algorithm  was  outlined  in  Section  3.2, 
Reference  1.  The  requirement  of  an  analytical  formulation  is  essential  to  WEDCOM, 
but  use  of  this  formulation  has  shown  some  unusual  reflection  altitudes  that  may 
not  be  appropriate.  Further  study  is  in  progress. 

Reflection  Coefficient  Computation  Criteria 
for  Adjacent  Vertical  Ionospheric  Profile 

Section  3.3,  Reference  1,  described  a  criterion  for  determining  reflection 
coefficient  sensitivity  to  changes  in  the  earth's  magnetic  field.  For  the  vast 
majority  of  executions  this  criterion  has  worked  well  in  minimizing  the  number  of 
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reflection  coefficient  computations.  However,  for  some  lightly  disturbed  night¬ 
time  cases  it  has  overestimated  reflection  coefficient  sensitivities  and  led  to 
too  many  computations.  This  algorithm  is  presently  under  study  for  the  final 
WEDCOM  V  Implementation. 

Reflection  Coefficient  Interpolation  Algorithm 

The  analytic  reflection  coefficient  representation  described  in  Section  3.4, 
Reference  1,  estimates  the  reflection  coefficient  using 

r^  *  exp  (A  +  BC^}  (3) 

where  is  the  cosine  of  the  ionospheric  angle  of  incidence  and  A  and  B 

are  the  complex  constants  computed  in  routine  BHRHS,  Figure  3.  This  procedure  has 
proven  very  reliable  for  interpolating  on  ionospheric  Incident  angle  for  most 
environmental  simulations.  It  unexpectedly  performed  very  well  for  anisotropic 
reflection  and  as  well  as  expected  for  isotropic  conditions.  This  formulation 
resulted  in  poor  estimates  of  reflection  coefficients  for  frequencies  greater  than 
200  kHz  and  for  VLF  cases  when  ref lection/absorption  occurred  over  a  thick  alti¬ 
tude  region.  The  first  problem  area  is  not  serious  since  most  cases  of  interest 
are  below  200  kHz.  The  second  problem  area  is  more  serious,  resulting  in  inac¬ 
curate  eigenangles.  The  problem  is  that  Equation  3  is  not  adequate  for  the 
simplified  eigenangle  algorithm.  The  solution  to  this  difficulty  is  described 
later. 

VLF  MODEL  FORMULATION 

The  reflection  coefficient  computation  criteria  for  adjacent  vertical  pro¬ 
files,  described  under  the  LF  model,  also  is  used  for  VLF,  with  appropriate 
frequency  dependence.  The  reflection  coefficient  interpolation  algorithm  also 
applies  to  both  the  LF  and  VLF  models,  but  with  some  special  modifications  for 
VLF.  The  following  subsections  summarize  submodel  development  which  applies  only 
to  the  VLF  model. 

Procedure  to  Reduce  the  Number  of  VLF 
Reflection  Coefficient  Computations 

WEDCOM  employs  a  modified  version  of  MODESRCH,  an  eigenangle  determination 
algorithm  developed  by  HOSC  and  described  in  Reference  8.  Section  4.2,  Reference 
1,  described  how  the  required  number  of  ionospheric  reflection  coefficient 


computations  could  be  reduced  by  25  to  50  percent  compared  to  WEDCOM  IV.  A 
further  reduction  In  the  number  of  ionospheric  reflection  coefficients  was  ob¬ 
tained  by  first  computing  reflection  coefficients  using  Equation  3  rather  than  for 
each  corner  of  the  complex-plane  area  being  searched.  The  reflection  coefficients 
are  computed  at  complex  ionospheric  incident  angles  with  real  part  defined  by 
ReCe^)  *  83,  81,  79,  77,  75  degrees,  and  imaginary  part,  Im(6i)  ,  that 
results  in  one  dB  attenuation  for  the  corresponding  value  of  Re(0^)  .  The  real 
and  imaginary  pairs  are  illustrated  in  Figure  5.  The  remainder  of  the  modified 
MODESRCH  procedure  remains  identical  to  that  described  for  WEDCOM  IV  in  Refer¬ 
ence  2. 

Eigenvalue  Determination  Using  Analytic  Reflection 
Coefficient  Formulation  and  Initial  Eigenangle  Estimates 

The  Initial  eigenangle  approximation  algorithms  described  in  Section  4, 
Reference  1,  were  Implemented  and  tested  for  the  matrix  of  cases  shown  in  Table  1. 
The  280-degree  azimuth  and  30-degree  dip  angle  correspond  to  strong  anisotropic 
conditions.  As  mentioned  earlier,  the  algorithm  worked  very  well  for  the  majority 


Figure  5. 


Illustrating  the  location  of  e.  used  for  reflection 
coefficient  computations. 


Table  1.  Parameters  varied  to  test  analytic  eigenangle  determinations. 


Frequency 

10,  15,  20,  25,  30  kHz 

Azimuth 

280  degrees 

Dip 

30,  70  degrees 

Profiles  = 

Normal  day  and  night 

spread  debris  parameter 

4  x  10-(20,  11,  10,  9,  8,  7,  6,  5,  4) 

Ground  Conductivity  = 

0.001,  0.01,  5.0  mho 

of  cases  including  highly  anisotropic  profiles.  However,  where  the  analytic 
reflection  coefficient  formulation  (Equation  3)  was  not  adequate,  a  few  important 
eigenangles  were  not  correctly  computed.  Several  procedures  were  attempted  to 
improve  the  eigenvalue  approximations.  The  improvements  were  partially  successful 
and  show  promise  of  complete  success.  If  time  permits  these  procedures  will  be 
pursued  on  the  present  effort  as  the  computational  time  reduction  and  simpler 
coding  makes  this  analytic  formulation  more  desirable  for  WEDCOM  than  the  search 
technique.  If  the  approximation  can't  be  completed,  a  brute  force  application  of 
the  MODESRCH  technique  will  be  Implemented  for  these  questionable  cases. 

Mode  Conversion  vs  WKB 

The  use  of  mode  conversion  provides  more  accurate  answers  than  the  WKB 
approximation  (an  averaging  technique),  when  propagation  conditions  vary  in  the 
direction  of  propagation.  It  "also  eliminates  the  need  for  defining  an  unambiguous 
mode  numbering  system.  However,  mode  conversion  requires  more  computation  time 
than  WKB.  For  this  reason  several  near-uniform  ionospheric  profile  cases  were  run 
using  both  the  mode  conversion  and  WKB  procedures  for  computing  field  strengths  as 
a  function  of  range.  The  comparisons  indicated  that  the  computationally  faster 
WKB  procedure  could  be  utilized  for  propagation  paths  with  small  variation  in 
ground  conductivity  and  ionospheric  profile  in  the  direction  of  propagation.  The 
use  of  WKB  when  appropriate  has  been  implemented. 


SECTION  4 

ELF  PROPAGATION  MODEL 


A  new  ELF  propagation  model  for  the  WEDCOM  code  was  described  in  Reference  1. 
The  model  was  adapted  from  a  model  developed  by  Booker  (Reference  9)  and  is  five 
to  ten  times  faster  running  than  the  current  WEDCOM  model.  Except  for  certain 
nighttime  ambient  ionospheres  noted  in  Reference  1  the  new  model  appears  to  give 
acceptable  accuracy  for  cases  of  interest.  In  the  following,  model  changes  to 
improve  results  for  ambient  nighttime  ionospheres  and  procedures  used  to  implement 
the  new  model  into  the  WEDCOM  code  are  described. 

REFLECTION  FOR  AMBIENT  NIGHTTIME  CONDITIONS 

Tables  2  through  5  show  phase  velocity,  c/v,  and  attenuation  rate,  a  , 
obtained  with  the  Booker  and  NOSC  (full-wave  solution)  models  for  ambient  night¬ 
time  profiles  (References  9  and  10).  Tables  6  and  7  show  the  electron  density, 
positive  ion  density,  and  collision  frequencies  used  in  the  two  nighttime  models. 
As  noted  in  Reference  1  the  results  are  in  good  agreement  for  the  nighttime 
profiles  used  by  Booker  but  show  differences  in  a  of  as  much  as  a  factor  of  two 
for  the  NOSC  nighttime  profile.  The  problem  is  related  to  the  electron  density 
used  in  the  E-region.  The  NOSC  profile  differs  from  that  used  by  Booker  in  that 
the  E-region  electron  density  minimum  is  too  large  to  cause  reflection  of  the 
extraordinary  (X)  wave  in  the  Booker  propagation  model  and  the  X  wave  energy 
reaching  the  E  region  is  not  returned  (thus  increasing  a).  Full-wave  solutions 
performed  at  NOSC  (Reference  10)  and  at  Tempo  show  that  considerable  energy  is 
returned  from  the  region  between  95  and  135  km.  The  calculation  also  shows  strong 
coupling  between  the  ordinary  (0)  wave  and  the  X  wave  in  the  E-region.  In  order 
to  include  the  energy  reflected  from  the  E-region  for  nighttime  conditions  the 
following  additional  computational  steps  have  been  added  to  those  described  in 
Reference  1. 
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Table  2.  c/v  and  a  obtained  with  Booker  model  for  the  Booker 
quiet  night  profile. 


Table  4.  c/v  and  a  obtained  with  NOSC  program  for  the  Booker 
quiet  night  profile. 


Table  5.  c/v  and  a  obtained  with  NOSC  program  for  the  NOSC 
quiet  night  profile. 


Table  6.  Profiles  used  for  comparisons  with  NOSC  calculations 


Alt  (km) 

Ne  (cm'3) 

N+  (cm-3) 

ve  (s) 

Vi  (s) 

0.00 

1  .99E-06 

5.501*03 

4 . 301*11 

1.076*10 

5.00 

3. 351-06 

4 . Qo£  *03 

1 .  JOE*  1 1 

5 . 30E  *09 

10.00 

6. 8JE-06 

3.50L+03 

8 .506*10 

2.606*09 

15.00 

l .551-05 

3.306*03 

3. 776*10 

l .296*09 

20.00 

3.546-05 

3.401*03 

1 . 6bL  *  10 

b. 376*08 

25.00 

8 . 45c— 0  5 

3 . 7C£ ♦ 03 

7. 4  5  6  *  J  9 

3.156*08 

30.00 

1.956-04 

3. 30£  *03 

3 . 306*09 

1.55E+08 

35.00 

4.006-04 

3 .OOE  *03 

1 . 476*09 

7 . 70E*  07 

40 .  00 

7 .676—04 

2.501*03 

b  •  5  0  E  *  0  6 

3.806*07 

4  5.00 

2.056— 0  3 

2.00E*03 

2 . 706 *oa 

1 .90E*u7 

50.00 

4.9«£-03 

l.50E*03 

l . 30E  *08 

9 . 30E  *  Ob 

55.  UO 

1.231-02 

l. 006*03 

5 . 706*07 

4 . b0£*06 

bO.OO 

3.211-02 

«. 006*02 

2 . but  *07 

2. 3QE+0b 

65.00 

5 . 57c— 02 

5 . OOE  *02 

1 . 13C*u7 

1.106*06 

70.00 

8.95E-02 

4 . 00  £ ♦ 02 

5 • OOE  *  06 

5.506*05 

75.00 

5 . OGE  +  0  0 

4 .OOE  *  02 

2.20E*06 

2.706*05 

ao.  oo 

7. 50E  *01 

3. 706*02 

9.96E+05 

1 . 34E ♦ 05 

85.00 

1 .606*02 

4 . OOE  *02 

4. 40E ♦ 05 

b.t4E  *04 

90.00 

3 . 30E  *02 

4. 50E*02 

l • 97E  *05 

3.236*04 

95.00 

6 .  JOE  -*0  2 

6 . OOE ♦ 02 

8. 701*04 

1.60E  *04 

10  0.00 

2.006*03 

2 • OOE  *  03 

3 . 90E  *04 

8.006*03 

105.00 

l.aOE*03 

1 • 80E  *03 

2 . bOE ♦ 04 

4 . Q0E*03 

110.00 

1 . 4GE*03 

1 . 40E  *03 

1 .806*04 

2 . OOE  *03 

115.00 

8.00E+02 

a. OOE *02 

1 . 04E  *04 

1.006*03 

120.0a 

6 . OOE  *  0 2 

6 . 00E+02 

1 .006*04 

3 . OOE  *02 

130.00 

3  •  706  *0  2 

3 . 70E ♦ 02 

4 . OOE  *0  3 

1 . 50E  *02 

14.0. 00 

2.30E+02 

2.806+02 

2 • 50E  *03 

8 . 006*01 

160.00 

2. 706*02 

2 . 70E  *02 

6.006*02 

1.806*01 

lao.oo  . 

4 . 50E  *0  2 

4 . 50E  *02 

1. 70E*02 

6.006*00 

200.00 

1 .0OE*03 

1 .00E*03 

4.506*01 

2.00E+00 

Table  7.  Booker  nighttime  profiles  (60-degree  latitude) 


Alt  (km) 

Nfi  (cm-3) 

N+  (cm-3) 

ve  (s) 

0.00 

l. 006-0 8 

3.156*03 

2.196*12 

10.00 

1 • 00c— 0  8 

3.196*03 

3.636*11 

20.00 

1 . 006-0 8 

3.096*03 

6 . 176*10 

30.00 

1.006-08 

2.756*03 

1 .  Oa  t  *  1 0 

90.00 

l. 006-08 

1. 9ib*0J 

1 . 786*09 

95.00 

l.GOfc-08 

1.906*03 

7.336*08 

50.00 

1.006-03 

9 . 9  3t*02 

3 . 02  6  ♦  0  8 

55.00 

l.OCc-08 

7 .056*02 

1.296*08 

60.00 

l.OOE-OS 

5.116*02 

5 .136*07 

65.00 

1 . COt— 0  8 

3.936*02 

2 . 116*07 

70.00 

1.26E-05 

3.176*02 

8 . 716*06 

75.00 

6.526-02 

2.976*02 

3.596*06 

80.00 

6.516*01 

2. 216*02 

1 . 95c *06 

85.00 

3  •  72t  *0  2 

9.376*02 

5.106*03 

86.00 

3.916*02 

9 . 92E* 02 

5 .116*05 

87.00 

9 • 06E  *0  2 

9.936-02 

9 . 2oc  *05 

88.00 

9 . 96c  *0  2 

9. 756*02 

3 .536  *05 

39. OC 

6.67  6*02 

6.896  *02 

3.006*05 

90.00 

1 . 92€-*0  3 

1 . 9  9fc ♦ 0  3 

2 .5lE*U5 

91.00 

2.126*03 

2.1 3fc*03 

2 .106*03 

92.00 

2.316*03 

2. 326*03 

1 . 7fc£ *03 

93.00 

2.356*03 

2 .36E*03 

l • 936*05 

99.00 

2.376*03 

2.3  7£*03 

1.296*05 

95.00 

2.386*03 

2.396*03 

L. 096*05 

100.00 

2.936*03 

2.936*03 

9.276*09 

105.00 

2.396*03 

2.396*03 

1 . 765*09 

110.00 

1.996*03 

1.996*03 

7.296  *0  3 

115.00 

1.196*03 

1. 1 96*03 

2.996*03 

120. OC 

5.016*02 

8.016*02 

1 • 23  £ *03 

122.50 

3.3  ofc  *0  2 

3 . 3  at  *02 

7 .  vuc.  *02 

125.00 

2 .36^*02 

2.386*02 

5.076*02 

127.50 

1.326*02 

1 . o  Zi *0  2 

j . 2pc  *02 

130.00 

1.516*02 

1.516*02 

2 . C96  *02 

190. 00 

1.396*02 

1.396*02 

3  .  t  *01 

195.00 

1 . 5 2  c *0  2 

1. 526*02 

1 . 9bc ♦ 01 

150.00 

1 . 806*02 

1. o  0  6  *  C i 

6 . 0 Jr ♦ Uu 

155.00 

2  .  lot  *  02 

2.1  6*  *  u. 

2 .  9s  l  ♦  0  0 

160.00 

2  .u  16  *3  2 

2 •  6  ic*  ».  J 

4  .  .  1  7U 

170.00 

3 . 8  3t*J  2 

3.836*02 

1  •  7  9  6  —  0  i 

180.00 

5.616*02 

■}.  tic  *02 

2 .956-02 

190.00 

3.296*32 

3.296*02 

5 .01:  —  03 

200.00 

1.216*03 

1.211*03 

5 .816—09 

vi  (s) 

3 . 39E* 1 1 
5.756*10 
9. 776*09 
1.666*09 
2  .82E*Oa 
1.  166*08 
9.796*0? 
1.976*07 
8. 136* 06 
3.35E*06 
1.38E*G6 
5.696*05 
2.396*05 
9 • 66 E  *09 
8 • 09  6  *  0  9 
6.7fc£  *09 
5.686*09 
9. 756*09 
3 . 98E  *09 
3.336*09 
2  •  79  E  ♦  0  9 
2.396*09 
1 . 96E*  0  9 
1.696*09 
6. 766*03 
2.796*03 
1. 15E*03 
9.736*02 
1.956*02 

I. 256*02 
6 . 09  £*01 
5.  lot  1 0 1 

J.  316*01 
5  .  6 2 6*00 
2 . 326*0 C 
9. 55c-0  1 
3.9ik-ol 
i . m2  t  — u  1 
2  .  756-0  2 
1.096-03 
7 . 99  £-0  9 
1  .  356-U9 
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Determine  the  E-region  altitudes  where  the  real  part  of  the  extraor¬ 
dinary  wave  index  of  refraction  reaches  maximum  and  minimum  values  (near 
100  and  140  km  for  the  NOSC  nighttime  profiles) . 

2.  Compute  the  equivalent  height  of  reflection,  Z^  ,  for  the  extra¬ 
ordinary  wave  for  reflection  heights  occurring  below  100  km  and  for 
E-region  reflection  heights  determined  in  Step  1.  The  phase  change  up 
to  the  bottom  E-region  reflection  height  is  taken  as  the  phase  change  up 
to  the  bottom  altitude  plus  0.63  percent  of  the  phase  change  between  the 
E-region  altitudes  determined  in  Step  1. 

3.  For  cases  where  the  Booker  model  does  not  result  in  reflections  above 

100  km  the  value  of  determined  in  Step  2  is  used  to  compute  the 

phase  velocity  and  attenuation  rate.  When  the  Booker  model  results  in 
reflection  above  100  km  the  attenuation  rate  is  computed  by  the  method 
defined  in  Step  2  and  by  methods  as  defined  in  the  original  model  and 
the  minimum  attenuation  rate  and  associated  phase  velocity  are  used. 

Tables  8  and  9  show  results  obtained  with  the  modified  model  for  the  two 
nighttime  profiles.  The  modification  does  improve  the  results  for  the  NOSC 
profile  and  either  does  not  change  or  slightly  improves  the  results  for  the  Booker 
profile. 


GROUND  LOSSES 

In  the  Booker  model  the  earth  is  taken  as  a  perfect  conductor  and  grouri 
losses  are  neglected.  Ground  losses  can  be  added  linearly  to  ionospheric  losses 
as  described  by  Wait  (Reference  11) 


a 


1  fQ  A  1  ]  .  ■ 

^  +  —  nepers  km 

2HR  N  Ng 


where 

a  =  attenuation  rate 

^  *  ionospheric  loss  term 
N 

Ng  =  ground  complex  index  of  refraction 
H„  =  effective  reflection  height  (km). 


Table  8.  c/v  and  a  obtained  with  revised  Booker  model  for 
the  Booker  quiet  night  profile. 


Latitude  (deg) 

Kiwis 

30 

45 

60 

75 

85 

75 

1.14 

1.13 

1.13 

1.13 

1.13 

1.13 

1.13 

0.39 

0.43 

0.44 

0.44 

0.44 

0.44 

0.44 

300 

1.08 

1.10 

1.11 

1.11 

1.10 

1.10 

1.10 

1.63 

3.37 

3.84 

3.98 

3.98 

3.77 

3.69 

Table  9 

.  c/v  and  a  obtained 

with  revised  Booker  model  for 

the  N0SC  quiet  night  profile. 

Latitude  (deg)  f 

mxSSmm 

bosks 

30 

45 

60 

75 

85 

75 

1.18 

1.22 

1.21 

1.20 

1.20 

1.19 

1.19 

1.38 

0.83 

0.86 

0.94 

1.02 

1.07 

1.09 

300 

1.11 

1.15 

1.11 

1.15 

1.19 

1.18 

mm 

3.11 

3.29 

3.86 

5.63 

3.87 

3.50 

In  terms  of  dB  per  Mm  the  ground  loss  term,  dB  ,  can  be  approximated  by 

„Bg  ;  8-72;  ‘O3  «!«-*  y- .■!  dB/Mm 

R 

where 

f  =  propagation  frequency  (Hz) 
a  =  ground  conductivity  (mho  m  ^). 

IMPLEMENTATION  OF  BOOKER  MODEL  INTO  THE  WEDCOM  CODE 
Computer  Routines 

The  following  routines  are  used  to  implement  the  Booker  model  into  the  WEDCOM 

code. 


Routine  ACCUMB 

This  routine  integrates  a  complex  function  f (x)  for  a  specified  interval  in 
x  assuming  either  a  linear  or  exponential  variation  with  x  .  Inputs  to  the 
routine  are  the  values  of  the  function  at  two  values  of  x  .  The  integral  of  the 
real  and  imaginary  parts  of  the  function  are  performed  separately.  If  the  value 


of  the  real  or  imaginary  part  of  the  function  changes  by  a  factor  of  2  or  less  a 
linear  variation  of  the  function  is  assumed;  otherwise  an  exponential  variation  is 
assumed . 

Routine  BFZ 

This  routine  computes  a  function  used  in  the  Booker  interpolation  scheme  to 
provide  smoothing  at  transition  altitudes  (Equation  48  in  Reference  1). 

Routine  BKRIAC 

This  routine  interpolates  a  function  of  the  complex  index  of  refraction  using 
a  four-point  interpolation  procedure.  The  derivative  of  the  function  can  also  be 
determined. 

The  four-point  interpolation  procedure  based  on  the  multiregion  fitting 

2 

technique  of  Booker  is  used  to  interpolate  the  function  n  -1  (see  Equations  48 
through  51  in  Reference  1).  Inputs  are  then  in  ( | A | )  and  in  ( | B | )  at  four 
altitudes  where  A  and  B  are  the  real  and  imaginary  parts  of  the  function.  The 
signs  of  the  real  and  imaginary  parts  of  the  interpolated  value  are  negative  for 
the  ordinary  wave  and  positive  and  negative  respectively  for  the  extraordinary 
wave.  The  wave  type  is  determined  from  the  input  MWAVE  (0  *  ordinary,  1  = 
extraordinary).  If  the  derivative  of  the  function  is  requested  (input  IDER  >  0), 
it  is  obtained  by  evaluating  Equation  51  in  Reference  1. 

Routine  BKRIAD 

This  routine  interpolates  a  real  function  using  a  four-point  interpolation 
procedure.  The  derivative  of  the  function  can  also  be  determined. 

Inputs  to  the  routine  are  either  the  function  or  the  log  of  the  function  at 
four  points  and  the  output  is  either  the  interpolated  value  or  the  log  of  the 
interpolated  value.  The  four-point  interpolation  procedure  based  on  the  multire¬ 
gion  fitting  technique  of  Booker  is  used  to  interpolate  the  function  (see  Equa¬ 
tions  48  through  51  in  Reference  1).  If  the  derivative  of  the  function  is  re¬ 
quested  (input  IDER  >  0)  it  is  obtained  by  evaluating  Equation  51  in  Reference  1. 

Routine  BOKELF 

This  routine  computes  the  complex  propagation  constant  at  a  specified  lo¬ 
cation  using  approximations  developed  by  Booker.  A  simplified  flow  chart  for  the 
routine  is  shown  in  Figure  6. 


First,  the  complex  index  of  refraction  (n)  is  computed  at  each  altitude  and 
2 

the  results  (actually  n  -1)  stored  by  calling  routine  RIQT.  For  daytime  con¬ 
ditions  calculations  are  only  made  at  and  below  100  km.  Next,  a  wave  loop  is 
established  to  perform  calculations  for  the  ordinary  wave  (MODEP=0)  and  the 
extraordinary  wave  (M0DEP=1).  For  each  wave  the  index  of  refraction  is  determined 
at  each  altitude  by  calling  routine  BKRIAC.  This  routine  determines  the  index  of 
refraction  using  a  four-point  fit  to  the  previously  stored  values  that  insures 
continuous  derivatives.  The  index  of  refraction  quantities  are  used  to  determine 
the  ratio  of  the  local  wavelength  to  the  refractive  index  local  scale  height.  For 
altitudes  above  95  km  the  altitudes  where  the  real  part  of  the  index  of  refraction 
for  the  extraordinary  wave  reaches  maximum  and  minimum  values  are  determined  (used 
in  computing  E-region  reflection  of  the  extraordinary  wave).  Then,  a  test  is  made 
to  determine  whether  there  is  a  reflection  altitude  in  the  last  altitude  interval. 
The  test  depends  on  the  reflection  altitude  number  and  the  wave  type  as  described 
in  Reference  1. 

If  there  is  a  reflection  altitude  within  the  interval  the  reflection  altitude 
is  computed.  An  iterative  procedure  is  used  to  account  for  the  change  in  re¬ 
flection  criteria  with  index  of  refraction.  After  the  reflection  altitude  is 
determined,  routine  BKRIAD  is  called  to  determine  ionization  and  collision  fre¬ 
quency  values  at  the  reflection  altitude  and  routine  RIQT  is  called  to  determine 
the  complex  index  of  refraction.  The  polarization  ratios  are  computed  (Equation 
36,  Reference  1)  and  the  phase  change  between  reflection  heights  determined 
(Equations  40  through  43,  Reference  1).  For  the  first  reflection  height  the 
equivalent  height  of  the  bottom  of  the  ionosphere  is  computed  (Equation  44, 
Reference  1).  In  computing  the  bottom  altitude  for  the  ordinary  wave  tests  are 
made  to  see  if  the  real  part  of  the  index  of  refraction  goes  to  zero.  If  it  does, 
the  height  where  it  goes  to  zero  is  determined  and  used  in  computing  the  bottom 
altitude. 

After  the  reflection  heights  and  other  quantities  have  been  computed  for  both 
the  ordinary  and  extraordinary  waves  the  equivalent  height  of  reflection  is 
computed  (Equation  39,  Reference  1).  For  nighttime  conditions  the  additional 
steps  described  in  this  report  are  used  to  account  for  nighttime  profiles  where 
the  extraordinary  wave  would  not  be  reflected  using  the  Booker  criteria.  Then  the 
complex  propagation  constant  is  determined  and  used  to  evaluate  the  phase  velocity 
and  attenuation  rate. 


The  mode  solution  is  determined  by  calling  routine  BOKELF  for  vertical  paths 
between  the  transmitter  and  receiver.  Then  ground  losses  are  determined  and  the 
mode  constant  adjusted  to  include  these  losses.  If  the  multiple  receiver  location 
option  is  exercised  (KALCMR  =  1),  excitation  and  height  gain  factors  are  found  for 
each  vertical  path  on  the  short  great  circle  path  between  transmitter  and  re¬ 
ceiver.  Otherwise,  excitation  and  height  gain  factors  are  only  found  for  the 
vertical  paths  at  the  transmitter  and  receiver  terminals.  If  field  strength 
calculations  for  propagation  along  both  the  short  and  long  great  circle  paths  are 
to  be  made,  ground  reflection  coefficients  and  excitation  factors  at  the  trans¬ 
mitter  and  receiver  terminals  are  computed  for  both  propagation  directions.  The 
excitation  and  height  gain  factors  for  the  long  great  circle  path  direction  are 
not  computed  for  receiver  locations  between  the  transmitter  and  receiver  terminals 
specified  in  input. 

Detailed  output  describing  the  eigenangle,  attenuation  rate,  relative  phase 
velocity,  ionosphere  reflection  coefficients,  ground  reflection  coefficients,  and 
excitation  factors  are  written  out  on  the  detailed  output  file  if  requested  (input 
option).  Finally,  routine  BOKTEM  is  called  to  evaluate  the  TEM  mode  field  values 
at  the  receiver. 

Routine  BOKTEM 

This  routine  calculates  TEM  mode  field  values  at  the  receiver  for  ELF  propa¬ 
gation.  A  simplified  flow  chart  for  the  routine  is  shown  in  Figure  7. 

First,  propagation  parameters  and  geometry  are  defined  and  nominal  system 
output  quantities  are  written  out.  Then  a  time  loop  over  the  calculation  times  is 
started  and  eigenvalue-dependent  data  parameters  defined  in  routine  ELFMOD  are 
read  from  a  data  file.  If  ambient  calculations  have  been  requested  (input  option) 
data  for  ambient  conditions  are  read  in  first;  otherwise,  data  for  disturbed 
conditions  are  obtained.  The  short  path  WKB  approximation  is  always  calculated 
and  the  long  path  calculation  is  made  when  requested  (input  option).  Then  the 
short  path  (and  the  long  path  when  requested)  electromagnetic  fields  at  the 
receiver  location  are  determined  and  nominal  output  is  written  out.  If  the 
multiple  receiver  location  option  (input  option)  has  been  exercised,  the  electro¬ 
magnetic  fields  at  each  vertical  path  location  are  computed.  After  calculations 
have  been  completed  for  disturbed  conditions,  calculations  are  started  for  the 
next  calculation  time. 


Figure  7.  Simplified  flow  chart  for  routine  BOKTEM 


Routine  ELFMOD 


This  routine  is  the  ELF  propagation  driver  routine  and  determines  anisotropic 
reflection  coefficients,  mode  solutions,  and  excitation  factors.  Only  the  trans¬ 
verse  electromagnetic  wave  is  analyzed.  A  simplified  flow  chart  for  the  routine 
is  shown  in  Figure  8. 

First  reflection  coefficients,  excitation  factors,  and  height  gain  parameters 
are  initialized.  Then  a  time  loop  over  the  calculation  times  is  started  and  an 
ambient  calculation  flag  (determines  whether  calculations  will  be  made  for  ambient 
conditions)  is  set  from  input  parameters. 

Next  a  loop  over  the  vertical  paths  between  transmitter  and  receiver  is 
started.  For  each  path  magnetic  field  and  ground  conductivity  quantities  are 
established  and  the  level  of  detailed  output  requested  determined.  The  ionization 
and  collision  frequencies  along  the  vertical  path  for  disturbed  conditions  are 
obtained  from  the  environment  file. 

For  ELF  propagation  there  are  NP-2  paths  on  the  great  circle  between 
transmitter  and  receiver.  When  calculations  are  performed  for  long  and  short 
paths,  only  NPT2  paths  are  along  the  short  path,  the  balance  along  the  long  paths. 
The  last  two  paths  (L  =  NP-1  and  L  =  NP)  are  located  a  Fresnel  zone  distance 
normal  to  the  short  great  circle  path  between  transmitter  and  receiver  and  are 
used  to  determine  whether  the  ionization  normal  to  the  great  circle  propagation 
path  is  spherically  stratified  (an  assumption  used  in  the  propagation 
calculations).  If  the  attenuation  rate  (dB  per  1000  km)  for  either  of  the  last 
paths  is  less  than  one  half  the  attenuation  rate  at  the  path  midpoint,  a 
diffraction  flag  is  set  to  indicate  that  the  propagation  calculations  may  be  in 
error. 


Routine  RIQT 

This  routine  determines  the  complex  index  of  refraction  at  a  point.  The 

index  of  refraction  (n)  is  computed  (Equation  31,  Reference  1)  and  the  quantity 
2 

n  -1  is  stored  for  the  ordinary  and  extraordinary  waves. 

Routine  SSQRT 

This  routine,  obtained  from  the  University  of  California,  San  Diego,  deter¬ 
mines  the  square  root  of  a  complex  number  with  a  consistent  sign  so  that  conti¬ 
nuity  in  defining  the  ordinary  and  extraordinary  waves  is  maintained  over  the 
entire  altitude  region  of  interest. 
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Figure  8.  Simplified  flow  chart  for  routine  ELFMOD 
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WEDCOM  Input  and  Output 

Input  to  the  WEDCOM  code  for  ELF  cases  is  the  same  with  the  new  model  as 
before  except  that  the  option  to  minimize  reflection  coefficient  calculations  is 
not  used.  Since  the  propagation  constant  calculations  are  much  faster  with  the 
new  model,  it  does  not  appear  efficient  to  eliminate  calculations. 

Most  of  the  output  from  the  ELF  model  is  identical  to  the  previous  model. 
Figures  9  and  10  show  examples  of  detailed  output  from  the  old  and  new  models 
respectively.  For  daytime  conditions  the  maximum  altitude  considered  is  100  km 
and  detailed  output  for  higher  altitudes  is  not  given.  In  the  old  model  the  index 
of  refraction  quantities  are  given  at  each  altitude  for  a  90-degree  magnetic  dip 
angle  (ordinary  wave).  The  new  model  results  show  output  for  both  the  ordinary 
(0)  and  extraordinary  (X)  waves  at  the  actual  magnetic  latitude  of  the  profile. 

An  additional  output,  ATTENG,  describing  the  ground  loss  in  dB  per  1000  km 
has  been  added  to  the  detailed  output  and  the  ground  reflection  coefficients 
previously  given  have  been  deleted.  Finally,  height-gain  factors  are  given  with 
the  excitation  factors  in  the  new  model  output. 
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Figure  9.  Illustration  of  ELF  detailed  output  for  old  model. 
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Figure  10.  Illustration  of  ELF  detailed  output  for  new  model. 


SECTION  5 

ATMOSPHERIC  NOISE  AND  SIGNAL  PROCESSING  MODELS 


This  section  describes  the  preliminary  work  required  to  implement  a  signal 
processing  capability  into  WEDCOM.  Previous  versions  of  WEDCOM  computed  signal 
propagation  parameters  such  as  the  received  EM  field  components  or  the  signal-to- 
noise  ratio  (S/N).  The  atmospheric  noise  magnitude  was  specified  by  the  user 
(using  CCIR,  Reference  12)  and  it  was  assumed  that  the  received  noise  value  was 
not  modified  by  a  nuclear  environment.  Estimates  of  message  quality  required 
additional  processing  not  provided  in  WEDCOM.  To  be  done  properly,  the  message 
quality  estimate  requires  knowledge  of  noise  characteristics  for  both  disturbed 
and  undisturbed  environments. 

This  additional  processing  required  was  awkward,  especially  for  analysts  with 
limited  propagation  and/or  signal  processing  background.  WEDCOM  V  will  incorpo¬ 
rate  a  signal  processing  model  which  includes  noise  sensitivity  to  nuclear  en¬ 
vironments. 

ATMOSPHERIC  NOISE  MODEL 

To  provide  useful  data  the  noise  model  must  define  the  characteristics  of  the 
received  noise.  Atmospheric  noise  for  the  lower  frequencies  is  characteristically 
a  homogeneous  background  superimposed  with  large,  intermittent,  nonoverlapping, 
pulses.  The  noise  is  due  to  thunderstorm  activity,  and  the  large  spikes  result 
from  local  lightning  strikes.  The  homogeneous  background  is  due  to  numerous 
unresolvable  pulses  from  distant  sources,  where  the  pulses  are  stretched  in  time 
due  to  propagation  dispersion. 

The  atmospheric  noise  can  be  characterized  as  received  noise  from  a  few 
transmitter  (noise  source)  locations.  This  assumption  is  especially  good  for 
receiver  modems  located  somewhat  distant  from  the  noise  sources.  For  example. 
Figure  11  illustrates  that  the  centroids  of  the  noise  sources  tend  to  be  centered 
south  of  the  equator  during  the  winter  season.  For  this  case  the  background  noise 
for  receivers  located  in  CONUS  could  be  modeled  by  a  single  noise  propagation 
path.  Winter  conditions  are  the  more  interesting  when  system  analysts  are 
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Figure  11.  Locations  of  winter  atmospheric  noise  radiation  centroids. 


concerned  with  maximum  noise  reduction  due  to  a  nuclear  environment.  For  summer 
conditions  (Figure  12)  the  noise  sources  tend  to  be  located  north  of  the  equator 
(closer  to  the  receiver  sites)  and  the  noise  modeling  now  requires  more  propaga¬ 
tion  links. 

Model  Requirements 

The  modem  models  (References  3  and  4)  require  that  the  two  following  noise 
characteristics  be  specified: 

1.  Average  noise  amplitude,  Vave  • 

2.  Ratio  of  the  noise  envelope  rms  to  its  average  amplitude,  . 

These  two  parameters  are  relatable  to  the  two  CCIR  parameters  (Reference  12).  A 
third  parameter  4*  ,  the  ratio  of  the  total  noise  power  to  the  background  Gaussian 
noise  power,  needs  to  be  specified.  Unfortunately,  V  is  not  a  CCIR  parameter 
and  is  difficult  to  predict.  Reference  3  provides  an  empirical  relationship  based 
on  experimental  data 

ndB)  =  2V,(dB)  U  -  [1.049/V,]2}  .  (4) 

d  d 

This  relationship  will  be  used  until  an  improved  prediction  becomes  available. 

To  satisfy  these  model  requirements  WEDCOM  V  will  have  at  least  two  of  the 
following  options: 

1.  Specification  of  the  receiver  noise  as  was  done  for  earlier  WEDCOM 
versions.  This  option  is  useful  for  the  cases  where  the  major  noise 
sources  are  local  or  are  not  due  to  thunderstorm  activity.  Minimum  code 
computations  are  required. 

2.  User  specification  of  the  noise  source  locations,  power  and  signal 
deviation. 

3.  Code  selection  of  the  noise  sources  based  on  user-specified  season  and 
diurnal  data.  This  option  requires  a  stored  data  base. 

The  ^ave  and  parameters  are  available  for  each  of  these  options.  The 

first  option  requires  the  CCIR  tables.  A  noise  prediction  model  is  required  for 
the  2nd  and  3rd  options,  to  be  executed  as  a  separate  activity  from  WEDCOM  or  as 
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Figure  12.  Locations  of  summer  atmospheric  noise  radiation  centroids. 
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Atmospheric  Noise  Prediction  Model 

An  atmospheric  noise  prediction  model  was  developed  by  the  Westinghouse 
Georesearch  Laboratory,  WGL  (Reference  13),  and  extended  to  include  elevated 
receiver  antennas  by  the  Naval  Research  Laboratory,  NRL  (Reference  14).  The 
WGL/NRL  code  was  developed  for  VLF  (10  to  30  kHz).  A  scaling  factor  developed  by 
the  Stanford  Research  Institute,  SRI  (Reference  15),  extends  the  frequency  range 
to  include  LF.  The  code  predictions  compare  well  with  extensive  worldwide  atmo¬ 
spheric  noise  measurements  for  vertical  polarization  and  with  limited  data  for 
horizontally  polarized  signals  received  by  elevated  antennas.  The  horizontal 
noise  component  is  related  to  anisotropic  conversion  from  vertical  to  horizontal 
polarization  at  the  receiver  location  only.  The  WEDCOM  propagation  models  are 
more  complete  than  the  NRL  models,  and  will  consider  quasi-TE  (as  well  as  quasi- 
TM)  modes  and  ionospheric  anisotropic  effects  along  the  entire  propagation  path. 
Thus,  the  horizontal  received  modes  (especially  important  to  A/C  with  trailing 
antennas)  will  receive  a  more  complete  evaluation.  A  limitation  to  the  WGL/NRL 
code  is  that  the  noise  source  is  assumed  to  be  vertical.  While  this  is  thought  to 
be  the  most  dominant  effect,  some  tests  (see  Reference  14)  have  shown  that  sig¬ 
nificant  horizontal  radiation  can  occur.  More  complete  modeling,  should  it  become 
available,  will  require  little  modification  to  WEDCOM  to  include  the  horizontal 
radiated  component. 

All  noise  sources  cannot  be  individually  evaluated  by  WEDCOM  and  an  algorithm 
which  will  minimize  the  number  is  required.  One  available  algorithm  is  found  in 
the  NRL  code.  The  noise  powers  are  combined  into  equivalent  noise  transmitters 
centered  in  5-degree  by  5-degree  grids  of  the  earth's  surface.  Where  possible, 
the  grids  are  further  combined  into  15-degree  by  15-degree  grids.  A  vector  sum  of 
the  propagation  modes  (or  rays)  is  not  meaningful  when  utilizing  this  algorithm  as 
the  noise  transmitter  locations  are  centroids  of  many  sources  and  the  timing  of 
all  sources  are  statistically  independent.  WEDCOM  will  therefore  only  compute  RMS 
values  of  the  received  noise. 

Nuclear  Environment  V _  and  Values 

ave  d 

The  prediction  of  ^ave  i-n  a  nuclear  environment  is  straightforward.  Once 
the  noise  powers  and  locations  have  been  specified,  propagation  links  to  the 
receiver  can  be  evaluated  and  received  noise  powers  determined. 

The  noise  impulsivity  parameter  is  not  easily  estimated  in  a  nuclear- 

disturbed  environment.  Two  techniques  have  been  used  by  analysts.  An  analytic 


formulation  is  provided  by  Field  and  Lewenstein,  Reference  16.  The  formulation  is 
based  on  a  single  link  uniform  propagation  path  for  the  Gaussian  background  noise, 
and  ignores  any  dispersive  propagation  changes.  Impulsivity  changes  occur  due  to 
suppression  or  enhancement  of  the  received  background  power. 

A  second  formulation  of  the  noise  statistics  was  presented  in  Reference  3  for 
use  with  the  modem  models.  For  mathematical  convenience  the  narrowband  atmo¬ 
spheric  noise  is  modeled  as  a  modulated  Gaussian  random  process,  rather  than  the 
sum  of  impulsive  and  background  noise  components  as  was  done  in  Reference  16.  The 
impulsivity,  ,  is  assumed  to  be  unmodified  by  the  nuclear  environment. 

Both  models  provide  good  comparisons  to  data,  but  a  direct  comparison  for 
WEDCOM  use  is  not  easily  determined,  and  an  extensive  effort  for  comparison  is 
not  warranted  due  to  the  lack  of  adequate  modeling  of  ¥  (Equation  4). 

Two  methods  of  predicting  V  in  a  nuclear  environment  are  presently  under 


consideration.  Each  method  is  a  variation  of  the  two  techniques  mentioned  above. 
In  the  first  method 


V,  is  estimated  in  the  following  manner: 
d 


1. 

2. 


Determine 

Determine 


ave 


V  and  ¥  for  the  normal  environment. 


(V  (disturbed)/!?  (normal)} 
ave  ave 


Ax  *=  { Vx (disturbed) /Vj (normal)  } 


where 
V 


3. 


1  is  the  received  signal  strength  from  the  nearest  noise  radiator. 
Determine: 

2 


-  (4) 


{¥  -  n  +  i 


4.  Determine  V,  in  disturbed  environment  by  solving  Equation  4  with  the 


substitution  of 


from  Step  3  above, 
a 


The  second  method  uses  a  variation  of  the  computation  that  is  part  of 

the  WGL  model: 

1.  The  N  most  significant  noise  sources  are  determined  in  the  undisturbed 

environment  (V  is  determined  from  these  N  sources).  The 

ave 

locations,  distances  from  the  receivers,  diurnal  modifier  index  and 
transmitter  powers  are  saved. 
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2.  From  the  N  sources,  the  propagation  paths  of  a  maximum  of  the  five 
most  significant  transmitters  are  evaluated  for  both  the  normal  and 
disturbed  environments. 

3.  is  determined  from  the  WGL  algorithm  (Reference  13),  but  with  the 
received  powers  from  Step  2,  rather  than  the  simplified  estimates  used 
in  the  WGL  model. 

Both  methods  will  be  evaluated  and  one  will  be  implemented  into  WEDCOM  V. 

SIGNAL  PROCESSING  MODEL 

The  Stanford  Research  Institute  (SRI)  developed  VLF/LF  modem  models 
(References  3  and  4)  for  the  MEECN  System  Engineering  Division  of  the  Command  and 
Control  Technical  Center  (CCTC) ,  Defense  Communication  Agency.  These  models 
estimate  modem  performance  in  an  atmospheric  noise  environment,  and  have 
demonstrated  a  specific  capability  for  predicting  performance  of  selected  MEECN 
modems.  The  model  allows  for  a  building  block  structure  of  simulating  non-linear 
noise  suppression  devices  (such  as  bandpass  clippers  and  hard  limiters)  combined 
with  various  types  of  generic  digital  demodulators. 

Figure  13  illustrates  a  simplified  receiver  block  diagram  that  incorporates 
the  individual  building  blocks  of  the  SRI  modem  model.  The  receiver  includes  one 
of  several  available  generic  matched-filter  models  with  or  without  a  bandpass 
clipper.  The  clipping  level  of  the  bandpass  clipper  must  be  user  specified.  The 
input  specifications  and  alternatives  are  listed  in  Table  10.  The  model  outputs 
include  the  bit-error-rate  (BER)  and  character-error-rate  (CER) ,  as  listed  in 
Table  11.  The  relationship  between  BER  and  CER  assumes  that  the  bit  errors  are 
independent.  For  an  N-bit  Baudot  encoded  character 

CER  =  1  -  (1  -  BER)N  .  (5) 

Error-correction  is  included  as  a  user  option,  i.e.,  a  15/11  Hamming  code  whose 
code  words  are  15  bits  in  length  and  contain  two  teletype  characters. 

A  working  version  of  these  models  has  been  implemented,  but  as  a  separate 
activity  from  WEDCOM.  Some  model  modifications  were  required  and  an  integration 
scheme  replaced.  However,  the  basic  model  remains  as  originally  implemented  by 
SRI.  Figure  14  illustrates  an  example  execution.  This  case  is  for  illustration 
and  does  not  pertain  to  an  existing  system.  The  effects  of  limiters  and  a  reduced 


42 


noise  environment  are  shown  (as  could  occur  in  a  nuclear  environment).  The 
*  1.049  dB  curve  corresponds  to  a  Gaussian  noise  background  where  the 
improvements  due  to  bandpass  clippers  are  minimal. 


Table  10.  Modem  model  input  specifications,  options  and  mnemonics. 


Definition 

Modem  Input 

Pseudo-noise  spread 

NO 

(none) 

spectrum  type 

CO 

(complementary) 

OR 

(orthogonal ) 

Modem  type 

PS  (phase-shift  key) 

MS  (minimal  shift  key) 

CF  (coherent  frequency  shift  key) 

CS  (compatible  shift  key) 

DC  (differentially  coherent  PSK) 

NC  (noncoherent  frequency  shift  key) 

Clipper  in  receiver? 

YE  (yes) 

NO  (no) 

Bit  rate  (per  second) 

0. <RATEsl000. 

Vd  bandwidth  (Hz) 

0.<BWVDsl0000. 

SNR  bandwidth  (Hz) 

0. <BNSRsl00000. 

Clipper  bandwidth  (Hz) 

Ratio  of  clipper  saturation 
level  in  signal-to-noise 

0. <BCs5000. 

rms  level  (dB) 

Processing  loss  due  to 

-60.sCSKs60. 

mismatched  filtering  (dB) 
Noise  impulsivity  parameter. 

O.sLOSSslO. 

Vd  (dB) 

Signal-to-noise  ratio  value 
(dB) 

1.049SVDS25. 

-60. sSNRs40. 

-60.sSNRs40. 


Table  11.  Modem  model  output  parameter  definitions. 

Efa/No  =  Bit-energy-to-noise-density  ratio 

BER  =  Bit-error-rate 

BCER  =  Character-error-rate  —  Baudot  code 

HCER  =  Character-error-rate  --  Hamming  15/11  code 

EFF  =  Clipper  signal-to-noise  ratio  enhancement  (dB) 
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Figure  13.  Block  diagram  of  simplified  receiver  model. 
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Interservice  Nuclear  Weapons  School 

DEPARTMENT  OF  THE  NAVY 

ATTN:  TTV 

Joint  Cruise  Missiles  Project  Ofc 

Joint  Chiefs  of  Staff 

ATTN 

JCMG-707 

ATTN:  C3S 

ATTN:  C3S,  Evaluation  Office  (HDOO) 

Naval  Electronic  Systems  Command 

ATTN 

PME  106-4,  S.  Kearney 

Joint  Data  Systems  Support  Ctr 

ATTN 

Code  501A 

ATTN 

G510,  G.  Jones 

ATTN 

PME  117-211,  B.  Kruger 

ATTN 

C-312,  R,  Mason 

ATTN 

Code  3101,  T.  Hughes 

ATTN 

C-500 

ATTN 

PME-106,  F.  Diederich 

ATTN 

G510,  P.  Bird 

ATTN 

PME- 1 17-2013,  G.  Burnhart 

ATTN 

PME  117-20 

National  Security  Agency 

ATTN 

W36,  0.  Bartlett 

Naval  Intelligence  Support  Ctr 

ATTN:  B-3,  F.  Leonard 

ATTN 

NISC-50 
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DEPARTMENT  OF  THE  NAVV  (Continued) 

DEPARTMENT  OF  THE  AIR  FORCE  (Continued 

Naval  Ocean  Systems  Center 

Rome  Air  Development  Center 

ATTN 

Code  532 

ATTN 

OCS,  V.  Coyne 

ATTN 

Code  5322,  M.  Paulson 

ATTN 

OCSA,  R.  Schneible 

ATTN 

Code  5323,  J.  Ferguson 

ATTN 

TSLD 

Naval  Research  Laboratory 

Rome  Air  Development  Center 

ATTN 

Code  4720,  J.  Davis 

ATTN 

EEP,  J.  Rasmussen 

ATTN 

Code  4700,  S.  Ossakow 

ATTN 

EEPS,  P.  Kossey 

ATTN 

Code  4108,  E.  Szuszewicz 

ATTN 

Code  7500,  B.  Wald 

Strategic  Air  Command 

ATTN 

Code  6700 

ATTN 

ADWA 

ATTN 

Code  4700 

ATTN 

DCX 

ATTN 

Code  4780 

ATTN 

XPQ 

ATTN 

Code  7950,  J.  Goodman 

ATTN 

XPFC 

ATTN 

Code  4187 

ATTN 

DCZ 

ATTN 

NRI/STINFO  Library 

Naval  Surface  Weapons  Center 

ATTN 

XPFS 

ATTN 

Code  F31 

OTHER  GOVERNMENT  AGENCIES 

Naval  Telecommunications  Command 

ATTN 

Code  341 

Department  of  Commerce 

National  Bureau  of  Standards 

Office  of  Naval  Research 

ATTN 

Sec  Ofc  for  R.  Moore 

ATTN 

Code  414,  G.  Joiner 

ATTN 

Code  412,  W.  Condell 

Department  of  Commerce 

National  Oceanic  4  Atmospheric  Admin 

Strategic  Systems  Project  Office 

ATTN 

R.  Grubb 

ATTN 

NSP-43,  Tech  Lib 

ATTN 

NSP-2141 

Institute  for  Telecommunications  Scient 

ATTN 

NSP-2722 

ATTN 

A.  Jean 

ATTN 

L.  Berry 

DEPARTMENT  OF  THE  AIR  FORCE 

ATTN 

W.  Utlaut 

Air  Force  Geophysics  Laboratory 

NATO 

ATTN 

LYD,  K.  Champion 

ATTN 

OPR-1 

NATO  Schoo 

(SHAPE) 

ATTN 

R.  Babcock 

ATTN 

US  Documents  Officer 

ATTN 

CA,  A.  Stair 

ATTN 

OPR,  H.  Gardiner 

DEPARTMENT 

OF  ENERGY  CONTRACTORS 

ATTN 

PHY,  J.  Buchau 

ATTN 

R.  O'Neil 

University 

of  Cal i fornia 

Lawrence  Livermore  National  Lab 

Air  Force  Weapons  Laboratory 

ATTN 

L-31,  R.  Hager 

ATTN 

SUL 

ATTN 

Tech  Info  Dept  Library 

ATTN 

NTN 

Los  Alamos 

National  Laboratory 

Air  University  Library 

ATTN 

T.  Kunkle,  ESS-5 

ATTN 

AUL-LSE 

ATTN 

D.  Simons 

ATTN 

MS  664,  J.  Zinn 

Assistant  Chief  of  Staff 

ATTN 

MS  670,  J.  Hopkins 

Studies  4  Analysis 

ATTN 

R.  Jeffries 

ATTN 

AF/SASC,  C.  Rightmeyer 

ATTN 

J,  Wolcott 

ATTN 

P.  Keaton 

Deputy  Chief  of  Staff 

Research,  Development,  4  Acq 

Sandia  National  Laboratories 

ATTN 

AFRDQI 

ATTN 

T.  Cook 

ATTN 

B.  Murphey 

Deputy  Chief  of  Staff 

Plans  and  Operations 

Sandia  National  Laboratories 

ATTN 

AFXOKCD 

ATTN 

Org  1250,  W.  Brown 

ATTN 

AFXOKT 

ATTN 

Org  4231  ,  T.  Wright 

ATTN 

AFXOKS 

ATTN 

Space  Project  Div 

ATTN 

D.  Thornbrough 

Electronic 

Systems  D1 v/SCT-2 

ATTN 

Tech  Lib  3141 

ATTN 

ESD/SCTE,  J.  Clark 

ATTN 

D.  Dahlgren 

Foreign  Technology  Division 
ATTN:  TQTO,  B.  Ballard 
ATTN:  NIIS  Library 


DEPARTMENT 

OF  DEFENSE  CONTRACTORS  (Continued) 

DEPARTMENT 

OF  DEFENSE  CONTRACTORS  (Continued) 

BDM  Corp 

Mission  Research  Corp 

ATTN 

L.  Jacobs 

ATTN 

R.  Bigoni 

ATTN 

T.  Neighbors 

ATTN 

G.  McCartor 

ATTN 

F.  Guigliano 

Berkeley  Rsch  Associates,  Inc 

ATTN 

F.  Fajen 

ATTN 

S.  Brecht 

ATTN 

D.  Knepp 

ATTN 

J.  Uorktnan 

ATTN 

Tech  Library 

ATTN 

C.  Prettie 

ATTN 

R.  Kilb 

ATTN 

R.  Bogusch 

University 

of  Cal  i  fornia  at 

San  Diego 

ATTN 

R.  Hendrick 

ATTN 

H.  Booker 

ATTN 

C.  Lauer 

ATTN 

S.  Gutsche 

Computer  Sciences  Corp 

ATTN 

R.  Dana 

ATTN 

F.  Eisenbarth 

Mitre  Corp 

EOS  Technologies,  Inc 

ATTN 

W.  Hall 

ATTN 

W.  Lelevier 

ATTN 

J.  Wheeler 

ATTN 

B.  Gabbard 

ATTN 

W.  Foster 

ATTN 

M.  Horrocks 

GTE  Communications  Products 

Corp 

ATTN 

R.  Steinhoff 

Pacific-Sierra  Research  Corp 

ATTN 

F.  Thomas 

GTE  Communications  Products 

Corp 

ATTN 

H.  Brode,  Chairman  SAGE 

ATTN 

J.  Concordia 

ATTN 

E.  Field,  Jr 

ATTN 

I.  Kohl  berg 

Physical  Research,  Inc 

Institite  for  Defense  Analyses 

ATTN 

R.  Deliberis 

ATTN 

H.  Gates 

ATTN 

J.  Devore 

ATTN 

J.  Aein 

ATTN 

T.  Stephens 

ATTN 

E.  Bauer 

ATTN 

J.  Thompson 

ATTN 

H.  Wolf hard 

R&D  Associates 

Johns  Hopkins  University 

ATTN 

R.  Turco 

ATTN 

C.  Meng 

ATTN 

M.  Gantsweg 

ATTN 

K.  Potocki 

ATTN 

C.  Greifinger 

ATTN 

J.  Phillips 

ATTN 

F.  Gilmore 

ATTN 

T.  Evans 

ATTN 

H.  Ory 

ATTN 

J.  Newland 

ATTN 

G.  Steyr 

ATTN 

P.  Komiske 

ATTN 

W.  Karzas 

ATTN 

W.  Wright 

Kaman  Sciences  Corp 

ATTN 

P.  Haas 

ATTN 

E.  Conrad 

Rockwell  International  Corp 

kaman  Tempo 

ATTN 

R.  Buckner 

ATTN 

B.  Gambill 

ATTN 

DAS I AC 

Rockwell  International  Corp 

ATTN 

W.  Schuleter 

ATTN 

S  .  Qu  flic  1 

ATTN 

U.  McNamara 

2  cy  ATTN 

W.  Knapp 

SRI  International 

2  cy  ATTN 

R.  Rutherford 

ATTN 

A.  Burns 

ATTN 

G.  Price 

Kaman  Tempo 

ATTN 

R.  Tsunoda 

ATTN 

DAS  I AC 

ATTN 

J.  Vickrey 

ATTN 

V.  Gonzales 

Maxim  Technologies,  Inc 

ATTN 

D.  Neilson 

ATTN 

E.  Tsui 

ATTN 

J.  Petrickes 

ATTN 

J.  Marshall 

ATTN 

M.  Baron 

ATTN 

R.  Morganstern 

ATTN 

R.  Livingston 

ATTN 

D.  McDaniels 

Mitre  Corp 

ATTN 

W.  Chesnut 

ATTN 

G.  Harding 

ATTN 

G.  Smith 

ATTN 

A.  Kymmel 

ATTN 

C.  Rino 

ATTN 

MS  J104/M,  R. 

Jresp 

ATTN 

W.  Jaye 

ATTN 

C.  Callahan 

ATTN 

R.  Leadabrand 

Science  Applications,  Inc 
ATTN:  M.  Cross 
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